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SUMMARY 

We  consider  the  equation 

ut  +  lxu  -  laf 1  (x)  P  °°udx  -  0,  u(x,Q)  ■  h(x), 

a  linearized  version  of  a  basic  aquation  in  plasma  physics. 
Deriving  a  relation  equivalent  to  the  usual  dispersion 
relation,  but  of  more  convenient  form,  we  derive  some  simple 
criteria  for  stability  of  solutions  of  this  equation. 
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OK  THE  STABILITY  OP  SOLUTIONS  OP  THE 
LINEARIZED  PLASMA  EQUATION 

Richard  Bellman 
John  M.  Richardson 


1.  Introduction 
The  equation 

(1)  ut  +  ixu  —  iaf '  ( x )  y7  00  udx,  u(x,0)  *  h(x), 

—00 

is  a  linearized  version  of  a  basic  equation  In  plasma  physics. 
Here  u  ■  u(x,t),  f(x)  is  a  real  function  of  x  and  a  is  a 
positive  constant.  A  problem  of  some  importance  is  that  of 
determining  when  the  solutions  of  (1)  are  uniformly  bounded  as 
t  — *•  oo  . 

The  usual  approach  to  this  problem  in  the  physical 
literature  is  to  write 

(2)  u-o8V(*). 

where  we  normalize  by  thr  condition  V700  /(x)dx  *  1.  Then 

c— oo 

(1)  yields 

(i) 

The  normalization  condition  yields  the  "characteristic 
equation"  for  8, 

<*> 


called  the  "dispersion  relation."  Since  it  is  generally  rather 
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difficult  to  study  the  distribution  of  the  roots  of  an 
equation  of  this  type,  we  wish  to  present  a  different  approach 
whloh  will  enable  us  to  present  some  simple  conditions  for 
stability. 

2.  Reduction  to  Renewal  Aquation 
Let  us  write  (1.1)  in  the  form 

(1)  ut  +  lxu  »  iaf '  (x)^00  udx,  u(x,Q)  -  h(x), 


and  regard  the  right-hand  side  as  a  forcing  term.  Then  we  may 
write 


(2) 


u  -  e 


-ixt 


h (x )  +  iaf  (x)  e”lx 8 ^ 

y0 


6/00u(x1,t)dx1 


—00 


ds  . 


Integrating  over  [-oo,ooj,  we  have 


(3)  /°°udx  -/°°e  lxth(x)dx 

-OO  ~ 


00  _ — lx  t  v 
-00 


/00f'(*)e-lx(t-a)dx 


-00 


p  00  udx. 


ds 


Let  us  now  Introduce  some  new  functions: 


(») 


v ( t )  -f°°  u(x,t)dx 
-00 


H( t )  - /°°  e  lxth(x)dx; 
—00 


K(t)  -y700  f '  (x)e  lxtdx. 


-00 


Then 
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(5)  v(t)  -  H(t )  +  /)t  K(t  -  a)v(s)de, 

"0 


an  integral  equation  of  renewal  type . 

The  Laplace  transform  readily  yields  the  equation 

<6> 


or ,  using  the  Laplace  inverse, 


where  C  can  be  taken  to  be  a  line  oarallel  to  the  imaginary 
axis  and  sufficiently  far  to  the  l  ight . 

The  question  of  the  asymptotic  behavior  of  u(x,t)  as 
t  — ♦  00  thus  reduces  to  the  problem  of  determining  the  root 
with  largest  real  part  of  the  equation 

(8)  1  -y’00 


3.  Discussion 

Using  the  definition  of  K(t)  given  in  (2.4),  we  see  that 
(2.5)  can  be  written 
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if  we  interchange  the  orders  of  integration.  However,  we 
definitely  do  not  wish  to  do  this.  The  advantage  of  using 
(2.5)  rather  than  the  dispersion  relation  lies  in  the  fact 
that  much  more  is  known  about  the  rootb  of  equations  of  the 
form  of  (2.3)  than  the  roots  of  equations  such  as  (1.4).  We 
shall  illustrate  this  in  the  next  section. 


4  .  Some  Simple  Criteria 

For  the  purpose  of  stability  investigations,  we  wish  to 
obtain  simple  conditions  which  ensure  that  (2.3)  possesses  no 
roots  with  nonnegative  real  part.  Although  a  simple  criterion 


a)  y700  |K(t)|dt  <  1, 
U0 


unless  K(t)  >  0  this  is  much  too  crude  for  useful  application. 

Let  us  suppose  that  K(t)  is  real,  and  write  s  ■  0  +  if*. 
Then  (2.3)  yields  two  equations, 

(2)  1  -  Z00  e"tft  cos  t?K(t)dt , 

% 

0  -  /  °°  e"K,t  sin  tTK(t)dt. 

°0 


Oiven  a  particular  function  K(t)  of  reasonable  analytic 
structure,  these  equations  can  be  used  to  determine  stability 
or  nonstability. 

A  particular  useful  and  simple  criterion  is  the  following 
Lemma.  If  g(t)  >0  and  monotone  decreasing,  then 
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{))  .P™  g(t)Bin  tmt  >  0 

% 

if  r  >0. 

Write 

(4)  J)  oo  g(t)8ln  tfidt  sin  ttTg(t )  -  g(t  ♦  J) 

0  ^0  L 

+  S(t  +  —  •  •  *Jdt . 

We  see  that  the  result  actually  holds  under  the  weaker  condition 

(5)  s(t)  -  g(t  +  I)  +  g(t  +  ^)  -  s(t  +  £)♦...>  0 
for  0  t  <  ¥/t . 

5*  An  Application 

As  an  application  of  the  foregoing  result,  let  us  consider 

_x2 

the  important  case  where  f(x)  ■  e  .  Since 

(1)  y  00  e-x2-lxtdx  .  ^  e-t2/\ 

—oo 

we  have 

(2)  y700  (-  21xe-x2)e-ixtdx  -  2V7  ^(e-t  /4) 

-  -  V7  te-t2/4 . 


Hence  the  equation  In  (2.5)  takes  the  form 

(3)  1  -  -  aVf,/00  te_ete_t2/ltdt. 

^0 

Integration  by  parts  yields 
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b 


(*) 


dt 


Ttius  the  equation  in  (})  becomes 


(5) 


1  +  2aVif 
2a^irs 


oo  -fit  -t/4 
e  e  ' 


dt . 


Write  s  ■  d  +  IV  and  equate  real  and  complex  parts .  The 
resulting  two  equations  are 


(6) 


1  +  2aVir 


2a  V/? 


JF7T2 

o  +  T 


■/ 


oo  _-dt 


e-t2/'> 


cos  trat, 


1  +  2aV? 


{  2a^¥ 


”'2"T72 

O  4-  t- 


00 


-K5tQ-t^/4 


sin  trdt 


Without  loss  of  generality  take  r  >  0.  Consider  first  the 

,  _ f-  ^  /h 

case  r  >  0.  Wien  if  d  >  0,  the  function  e  e  '  ie 
positive  and  monotone  decreasing  for  t  ^  0.  Hence  the  second 
equation  in  (6)  cannot  hold  for  a  nonnegative  value  of  d. 

If  t  -  0,  the  first  equation  yields 


(7) 


1  +  2aVi 


l  2a\/»  J 


w 


co 


-<tte— tV4 


dt. 


a  relation  which  holds  for  no  value  of  d  since 


(8) 


1  +  2a V? 


2e.Vt 


a-/ 


00  -dt  -tv1* 


dt  </ 


00  1 

e  dt  •  -r. 

d 


In  the  case  of  general  f(x),  we  obtain  criteria  in  terms 
of  the  Pourier  transform  of  f'(x),  or  of  f(x).  In  this  way 
we  can  exhibit  families  of  functions  which  yield  stability. 


